Helicoidal magnetic order in a clean copper oxide spin chain compound 
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We report susceptibility, specific heat, and neutron diffraction measurements on NaCu2 02, a 
spin-1/2 chain compound isostructural to LiCu202, which has been extensively investigated. Below 
12 K, we find a long-range ordered, incommensurate magnetic helix state with a propagation vector 
similar to that of LiCu2C>2. In contrast to the Li analogue, substitutional disorder is negligible in 
NaCu202. We can thus rule out that the helix is induced by impurities, as was claimed on the basis 
of prior work on LiCu2 02- A spin Hamiltonian with frustrated longer-range exchange interactions 
provides a good description of both the ordered state and the paramagnetic susceptibility. 

PACS numbers: 75.10.Pq, 75.40.Cx, 75.25.-fz 



Copper oxides are excellent model systems for low- 
dimensional spin-1/2 quantum antiferromagnets. In par- 
ticular, copper oxides with magnetic backbones com- 
prised of chains of CuO.4 squares have been shown to 
exhibit quasi-one-dimensional behavior. Two classes of 
copper oxide spin chain materials are known. Com- 
pounds in which adjacent squares share their corners are 
excellent realizations of the one-dimensional (ID) spin- 
1/2 Heisenberg Hamiltonian Q, [H 01 Linear Cu-O-Cu 
bonds along the spin chains give rise to a large antifcr- 
romagnetic nearest-neighbor exchange coupling. In com- 
pounds built up of edge-sharing squares, on the other 
hand, the Cu-O-Cu bond angle is nearly 90°, so that the 
nearest-neighbor coupling is more than an order of mag- 
nitude smaller jij. Because of the anomalously small 
nearest-neighbor coupling, longer-range frustrating ex- 
change interactions have a pronounced influence on the 
physical properties of these materials. Edge-sharing cop- 
per oxides thus provide uniquely simple model systems 
to test current theories of spin correlations in frustrated 
quantum magnets. 

At low temperatures, the ground state of edge-sharing 
copper oxides is either a 3D-ordered antiferromagnet 
[ElI or a spin-Peierls state || , depending on whether 
interchain exchange interactions or spin-phonon interac- 
tions are dominant. In the former case, the magnetic or- 
der is almost always collinear. An interestin g ex ception 
was recently discovered in LiCu 2 2 0, fiol Tllj . which 
undergoes a transition to a magnetic helix state at low 
temperatures. While such a state is expected for clas- 
sical spin models with frustrating interactions, quantum 
models predict a gapped spin liquid state in the range 
of exchange parameters that was claimed to describe the 
spin system in LiCu2 02- Since the ionic radii of Li + and 
Cu 2+ are similar, chemical disorder was identified as a 
possible solution to this puzzle. Indeed, a chemical anal- 
ysis of the sample used in the neutron scattering study 
of Ref. 0| showed that about 16% of the Cu 2+ ions in 
the spin chains were replaced by nonmagnetic Li + impu- 
rities. Since even much lower concentrations of nonmag- 
netic impurities are found to induce magnetic long-range 
order in other quasi- ID spin-gap systems, the authors 
of Ref. 0] attributed the unexpected helix state to the 



Here we report magnetic susceptibility, specific heat, 
and neutron diffraction data on NaCu2 02, a Mott insula- 
tor that is isostructural to LiCu2C>2. However, since Na + 
is much larger than Cu 2+ , substitutional disorder is a pri- 
ori unlikely in NaCu202- Chemical analysis and neutron 
diffraction data confirm that Na-Cu inter-substitution is 
negligible in our NaCu2C>2 samples. Below a Neel tem- 
perature Tn of 12K, we find an incommensurate helix 
state similar to that in LiCuaOq. Contrary to the inter- 
pretation advocated in Ref. [11J, this state is thus the 
ground state of a spin chain system without impurities. 
The incommensurability along the chain axis is such that 
the magnetic unit cell is nearly quadrupled with respect 
to the chemical cell. Together with an analysis of the sus- 
ceptibility data, this indicates that an antiferromagnetic 
next-nearest-neighbor interaction along the spin chain is 
the dominant exchange interaction in NaCu202- This 
disagrees with the spin Hamiltonian recently proposed 
for LiCu2C>2 0| . Each copper oxide chain thus con- 
tains two interpenetrating, but nearly independent ID 
spin systems. A model with longer-range exchange inter- 
actions provides a quantitative description of the ground 
state and paramagnetic susceptibility of this spin-1/2 sys- 
tem. 

Micro-crystalline powder samples of NaCu202 were 
synthesized via the azide /nitrate route in specially de- 
signed containers [l2l Il3 . fl4| . The starting materials 
were milled, pressed in pellets under 10 5 N, dried in vac- 
uum (10~ 3 mbar) at 150° C for 12 h, and placed in an 
argon atmosphere in a tightly closed steel container pro- 
vided with a silver inlay. Finally, in a flow of dry argon 
the following temperature sequence was applied: 25 to 
260° C (100 °C/h); 260 to 380° C (5 °C/h), 380 to 500°C 
(50 °C/h) with subsequent annealing for 30 hours at 500° 
C. Powder X-ray diffraction patterns for initial charac- 
terization were collected with a Stoe STADI-P diffrac- 
tometer using CuK Q i and MoK Q i radiation. No impurity 
phases were detected, except small traces of CuO and 
CU2O. The single crystals were grown by the self- flux 
method in an argon atmosphere in platinum crucibles 
(lif . The crystals have a platelet shape with typical sizes 
of 7 x 3 mm 2 and thickness of up to 100 [im. X-ray 
diffraction from single crystals ground to powder showed 
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FIG. 1: (a) Crystal structure, and (b) proposed magnetic 
structure of NaCu2 02. In panel (b), the most important mag- 
netic interactions within a chain pair are sketched. 



In addition, no platinum impurities were detected in our 
crystals by inductively-coupled plasma atomic emission 
spectroscopy (ICP-AES) measurements. Finally, we have 
used ICP-AES to establish the chemical composition of 
our crystals as Na1.ooio.02Cu2.ooio.02Oj,. The Na:Cu ra- 
tio is thus identical to the ideal stoichiomctry within the 
experimental error. This should be compared to the sam- 
ples studied in Rcf. 1 1] whose chemical composition was 
quoted as Li1.i6Cu1.s4O2. 01 ■ 

For the neutron diffraction measurements, a powder 
sample of weight ~ 5 g was fitted into a cylindrical 
vanadium container and loaded in a standard helium- 
flow cryostat. High-resolution diffraction patterns were 
collected at 3.5, 23 and 300 K on the two-circle pow- 
der diffractometer D2B at the Institut Laue-Langevin, 
Grenoble. A wavelength A = 1.5946A was selected by 
a Ge[355] monochromator. The diaphragms were set 
to 70/70 and the collimations to 10'-10'-10' before the 
monochromator, between monochromator and sample, 
and between sample and detector, respectively. Com- 
plementary neutron diffraction data were taken on the 
high-flux powder diffractometer D20 using a wavelength 
A = 2.4A. The diaphragms were set to 70/10. Mea- 
surements of the magnetic susceptibility and the heat 
capacity were carried out on powder samples as well as 
on single crystals. The powder measurements were taken 
on a sample of total mass 115 mg in a SQUID magne- 
tometer (MPMS, Quantum Design). 

A Rietveld analysis of the diffraction pattern at room 
temperature confirms that NaCu202 has an orthorhom- 
bic crystal structure (space group Pnma) with lattice pa- 
rameters a = 6.2087A, b = 2.9343A and c = 13.0548A, 
in agreement with an X-ray investigation of small single 
crystals. A pictorial representation of the lattice struc- 
ture is shown in Fig. 1. The unit cell of NaCu202 con- 
tains four magnetic Cu 2+ ions belonging to two pairs of 
copper oxide chains running along b. The Cu-O-Cu bond 
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lent angle in LiCu2 02- The chains are separated from 
each other in the a-direction by rows of non-magnetic 
Cu + ions, and in the c-direction by Na + ions. Since the 
two chains within one pair are shifted relative to each 
other by half a unit cell along b, each chain can also be 
viewd as a single zig-zag chain [lCI . The nearest- neighbor 
Cu 2+ -Cu 2+ distances between different chain pairs are 
considerably larger than those within a pair. It is thus 
reasonable to expect that the system is magnetically one- 
dimensional. 

In order to check for possible Na-Cu inter-substitution, 
we have refined the occupancies of Na and Cu sites in the 
Rietveld analysis of the neutron diffraction data and ob- 
tained full occupancy of both positions within an error of 
2%. This finding is expected based on the different ionic 
radii of Na + and Cu 2+ , and it confirms the results of the 
chemical analysis above. In contrast to LiCu2 02, substi- 
tutional disorder on the copper oxide chains is therefore 
negligible in NaCu2 02- 

Fig. 2 shows the uniform magnetic susceptibility of 
a NaCu2 02 powder as a function of temperature. The 
main feature of the curve is a broad maximum at 52 K, 
which is characteristic of low-dimensional spin systems 
and indicates a crossover to a state with antiferromag- 
netic short-range order. A Curie- Weiss fit of the high- 
temperature susceptibility for 200 < T < 300 K yields 
a negative Curie- Weiss temperature of @cw = —62 K, 
indicating predominant antifcrromagnetic interactions. 
Because of the intrinsic ID nature of our system, we have 
fitted the susceptibility curve to the exact solution of the 
5 = 1/2 Heisenberg chain with a single antiferromagnetic 
interaction J2 |l7j |. resulting from the dominant Cu-O- 
O-Cu exchange path. We added a diamagnetic contri- 
bution from the closed atomic shells, which we estimate 
as -52 xl0 _6 cm 3 /mol from Pascal's increments The 
fit yields J2 = 85 K (fcs=l) and g — 2.07, and the re- 
sult is shown in Fig. 2. Our susceptibility data are thus 
in good quantitative agreement (except small deviations) 
with a model including a single antiferromagnetic interac- 
tion parameter along the spin chains (Fig. 2), supporting 
the view of a single chain as two interpenetrating, weakly 
coupled J 2 -Heisenberg chains. 

The low-temperature susceptibility and specific heat 
data plotted in Fig. 3 indicate two magnetic phase tran- 
sitions at 12 and 8 K. Both transition temperatures are 
much smaller than J2, as expected based on the quasi- 
one-dimensionality of the magnetic lattice. The temper- 
ature dependent susceptibility of a single-crystal sam- 
ple in the three crystallographic directions is plotted in 
Fig. 3a. The susceptibility in the 6-direction, \b, is not 
strongly affected by the 12 K transition, while Xa and Xc 
are suppressed. Below the 8 K transition, on the other 
hand, Xa and Xb exhibit an upturn, reflecting the role 
of anisotropic interactions. The specific heat anomaly at 
the 8 K transition is weaker and obliterated by a modest 
magnetic field of 9T (Fig. 3b), whereas the 12 K transi- 
tion is more robust. Taken together, these data suggest 
that the 8 K transition arises from spin canting. A simi- 
lar rantino' transition has hppn nhsprvpd in TinCiiDn I7K 
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FIG. 2: Temperature dependence of the static magnetic sus- 
ceptibility of a powder sample of NaCu2 02 in a magnetic field 
B=0.1 T (open circles). The thick red line is the result of a 
fit to the helix model described in the text, whereas the blue 
line (dot-dashed, mostly hidden by the red one) is the (ex- 
act) susceptibility of the lD-Heisenberg chain (J2=85K) (l7| . 
For comparison, the green line represents the result of a fit 
to a simpler J1-J2 helix model (see text). The inset shows 
the actually measured data, from which a Curie contribution 
XCur=0.0105/T has been subtracted to obtain the data pre- 
sented in the main panel. Note that the Curie tail is very 
weak, providing testimony to the high quality of the sample. 




FIG. 3: Single-crystal data for NaCu202 as a function of 
temperature, (a) Susceptibility taken with a magnetic field 
of 0.1T applied along the principal crystallographic axes, (b) 
Specific heat for two applied magnetic fields: 0T (open cir- 
cles) and 9T (full circles) . The inset extends the temperature 
range. 



These data indicate that the magnetic ground state 
is more complex than suggested by the initial analysis 
of the paramagnetic susceptibility. In order to deter- 
mine the spin configuration in the ground state, neutron 
powder diffraction data were taken on the high-flux pow- 
der diffractometer D20. The results are shown in Fig. 
4. Several additional Bragg reflections are observed be- 
low 12 K at low scattering angles. Since no high-angle 
counterparts of these reflections are observed, and the 
intensities vanish at the magnetic ordering temperature 
(inset in Fig. 4), these reflections can be identified as 
magnetic. They can be indexed based on the magnetic 
propagation vector (0.5, £, 0) with ( = 0.227, which cor- 
responds to a pitch angle 4>b — 81.7° along the 6-axis. 
Since an amplitude modulation of the copper spin is not 
expected in a Mott insulating state, the intensities were 
compared to models in which the spins form a helicoidal 
state with identical amplitude on every Cu 2+ site, anal- 
ogous to that observed in LiCu2 02. An excellent refine- 
ment of the peak intensities (solid line in Fig. 4) was 
obtained by assuming a helix polarized in the a6-plane 
(Fig. 1), similar to LiCu 2 2 [Hi [13. The ordered mo- 
ment per Cu 2+ ion at base temperature is 0.56(4)/^. 

Since the observed helical spin structure of NaCu2 02 
cannot be accounted for within a simple nearest-neighbor 
Heisenberg model, we reanalyzed the paramagnetic sus- 
ceptibility in terms of a spin Hamiltonian with longer- 
ranged interactions 3d- Such additional couplings are 
expected for Cu-O-O-Cu bridges in edge-sharing copper- 
chains 0|. We consider magnetic couplings up to a dis- 
tance d=4b. Whereas couplings Jd>2 are always antifer- 
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magnetic depending on the Cu-O-Cu bond angle [4J. The 
classical ground state of this model is a helix whose pitch 
angle (f>t, is given by the expression 



cos (pb 



Ji - 3 J 3 
~4(J 2 -3J 4 ) ' 



(1) 



which is valid as long as J\ | , | J3 and | J4 1 are significantly 
smaller than \3i\. Using a finite temperature Lanczos 
technique |2jj, we have calculated the temperature de- 
pendent susceptibility on chains of length N — 24 for 
models where condition Q is fulfilled, using the angle 
4>b = 82° obtained in the analysis of the neutron scat- 
tering data. An excellent fit (red line in Fig. 2) is ob- 
tained for the entire temperature range from just above 
Tn up to 330 K, with the parameters J\ = —16.4 K, 
J 2 = 90 K, J 3 = 7.2 K, and J 4 = 6.3 K. The g-factor, 
g = 2.14, is typical for Cu 2+ in square-planar geome- 
try. The fact that the antiferromagnetic interaction J2 
is by far the largest parameter in the spin Hamiltonian 
explains the surprisingly good fit obtained by the nearest- 
neighbor model discussed above. The longer-range cou- 
plings J3, J4 are small, as expected, yet they are essential 
in this context, since fits attempted with J3 = J4 = 
(J1/J2 = —0.52 which yields the experimental pitch an- 
gle, green line in Fig. 2) are in significantly poorer agree- 
ment with the experimental data. 

We have thus determined the magnetic ground state 
of NaCu202, a spin-1/2 chain compound without signif- 
icant substitutional disorder, as a magnetic helix with a 
substantial ordered moment of 0.56/is. This rules out 
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FIG. 4: Difference of neutron diffraction data collected at 2K 
and 20K (open circles). The extra reflections at low temper- 
ature are all of magnetic origin. The solid line is a fit to a 
spiral model in the bc-plane with propagation vector (0.5, £, 0) 
with £ = 0.227. The temperature dependence of the (0.5, 1) 
reflection is plotted in the inset together with a fit of the in- 
tensity to a power law. 

proposed ^lj f° r LiCu202, where a helix with a similar 
propagation vector was observed. The spin Hamiltonian 
is dominated by a large antiferromagnetic exchange pa- 
rameter Ji ~ 90 K, which can be identified as the inter- 
action between second-nearest neighbors along the spin 
chains. The dominance of the J2 coupling may at first 
seem surprising, but it can be understood as a conse- 
quence of the structural geometry of NaCu202, whose 
Cu-O-Cu bond angle, 92.9°, is very close to the critical 
angle of ^94° at which the nearest-neighbor interaction is 
expected to change sign Q. Hence, the nearest- neighbor 
coupling Ji is small. The magnitude of J2 is in very good 
agreement with theoretical calculations for edge-sharing 
copper oxide chains with Cu-0 bond parameters deter- 
mined for NaCu2C>2 @]. This agreement is significant, 
because the prediction is largely insensitive to the bond 
angle and hence quite robust, and because direct mea- 
surements of this parameter in other copper oxide spin 
chains yield comparable results. The inclusion of a small 
interchain interaction within a chain pair (J 7 in Fig. lb) 
marginally affects the values of the intrachain exchange 
parameters extracted from the susceptibility, and there- 
fore does not change this picture qualitatively. 

In view of our work on a clean, isostructural compound 
and the theoretical work of Ref. pj], the exchange pa- 
rameters of LiCu2 02 should be reexamined. Because 
of the smaller Cu-O-Cu bond angle, the magnitude of 
the nearest-neighbor coupling in the spin Hamiltonian for 
LiCu2C>2 is expected to be larger than that of NaCu202, 
but the next-nearest-neighbor coupling should be com- 
parable. A scenario in which an antiferromagnetic in- 
terchain interaction of magnitude J' ~ 70 K significantly 
exceeds all of the relevant intrachain interactions 0] ap- 
pears very unlikely |22| . 



In conclusion, the magnetic helix state and paramag- 
netic susceptibility of NaCu2C>2 are in good quantita- 
tive agreement with the predictions of a model including 
longer-range exchange interactions. An open question 
concerns our use of the classical expression Q for the 
pitch angle. Some theoretical work suggests that the 
pitch angle for quantum models deviates substantially 
from this expression |2lj . However, one has to keep in 
mind that the magnitude of the ordered moment per 
magnetic copper site is 0.56/is, somewhat larger than 
the values typically observed in 1-D systems Quan- 
tum zero-point fluctuations thus appear to be signifi- 
cantly suppressed in the ordered state. Large ordered 
moments were also observed in other corner-sharing cop- 
per oxide spin chain compounds, and were ascribed to 
large exchange anisotropies. A quantitative understand- 
ing of the interplay between spin anisotropies and frus- 
tration in NaCu2C>2 and other spin chain materials is an 
interesting subject of further investigation. 
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